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Stochastic equation of fragmentation and branching 
processes related to avalanches 


Lucian Beznea 1 , Madalina Deaconu 2 , and Oana Lupa§cu 3 

Abstract. We give a stochastic model for the fragmentation phase of a snow avalanche. 
We construct a fragmentation-branching process related to the avalanches, on the set of all 
fragmentation sizes introduced by J. Bertoin. A fractal property of this process is emphasized. 
We also establish a specific stochastic equation of fragmentation. It turns out that specific 
branching Markov processes on finite configurations of particles with sizes bigger than a strictly 
positive threshold are convenient for describing the continuous time evolution of the number 
of the resulting fragments. The results are obtained by combining analytic and probabilistic 
potential theoretical tools. 
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1 Introduction 

The snow avalanches were studied from different points of view. Deterministic avalanche mod¬ 
els of Saint-Venant type, including numerical simulations, were investigated, e.g. [PuHu 07], 
[WaHuPu 04], [GrAn 09], [BouFe 08], [IoLu 15], and [IoLu 15a|. 

Fractal properties were also emphasized in a natural way, including probabilities. A fractal 
model for grain size distribution of a snow avalanche, by introducing the concept of aggregation 
probability, as a coagulation mechanism, is developed in [FaLoGr 03]. The fractal character 
of the snow has been studied in the paper |DeBiCh 12], In [CaChFr 10] are investigated the 
scaling properties of snow density in a stochastic fractal framework, able to reproduce the local 
randomness of real microstructure. 

There are several attempts to relay the rock avalanches to some fragmentation processes; cf. 
[DeBl 11], [DeBl 14], and [GrAn 09], however, the used fragmentation process is deterministic. 
In the articles [GaSaSc 93], [Lee 04], and [Za 95] a discrete time branching process is associated 
to an avalanche. 

In this paper we give a stochastic model for the fragmentation phase of a snow avalanche. 
The continuous time evolution of the number of the resulting fragments is described by using 
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specific branching Markov processes. More precisely, using a recent approach for the fragmen¬ 
tation processes developed in |BeDeLu 15], we construct a fragmentation-branching process 
related to the avalanches, on the set of all fragmentation sizes (:= the set of all decreasing 
numerical sequences bounded above from 1 and with limit 0; cf. [Ber 06]). A fractal property 
of this process is proved. We use analytic and probabilistic methods from the potential theory 
associated to semigroups and resolvents of kernels. 

The description of our avalanche model is the following. We fix a rupture factor r G (0,1), 
corresponding to an uniform proportionality of the fragments, constant in time. This specific 
property of an avalanche will be encoded in our given discontinuous fragmentation kernel. A 
Markov chain is induced in a natural way: the state space is E = [0,1], regarded as the space 
of all fragment sizes, from a point x G E it is possible to move to the point fix and (1 — /3)x 
with the respective probabilities (3 and 1 — /?, where fi := j 2 —. It is an immediate observation 
that if we assume fi < ~ 2 then the probability of occurring a bigger fragment is bigger then a 
smaller one. A time continuous Markov chain (actually, a continuous time random walk) is then 
canonically constructed, it is an analog of the Markov process constructed from a fragmentation 
equation associated to a continuous fragmentation kernel; see |BeDeLu 15] for details and other 
references. 

In the next section we recall (for the reader convenience) some properties on the compound 
Poisson pure jump processes and we introduce the discontinuous fragmentation kernels for 
avalanches we need further to develop our stochastic model. We fix a sequence of thresholds 
{d n ) n ^ i C (0,1) strictly decreasing to zero. In Proposition 2.2 we establish compatibility 
properties between the continuous time Markov chains X n , n ^ 1, constructed on each segment 
[. d n , 1]. In addition, every process X n enables to solve a corresponding martingale problem, cf. 
assertion ( ii ) of Remark 2.1. Note that we cannot associate a stochastic differential equation 
of fragmentation with the method from |FouGi 03] and |BeDeLu 15], because in the present 
context the fragmentation kernel is discontinuous. However, in Section 3 we write down in 
Theorem 3.1 the corresponding stochastic equation of fragmentation, that is, having the solution 
equal in distribution with X n . We use essentially the already mentioned existence of the 
solution of the martingale problem and the technique from [Jac 79] to derive a solution for 
the associated stochastic integro-differential equation, which will eventually lead to the claimed 
stochastic equation of fragmentation. 

As for the case of the fragmentation equation, we are interested in studying the evolution in 
time of the number of fragments of the avalanche, having the same size. In Section 4, to handle 
such a problem, we first investigate the time evolution of fragments bigger then a given strictly 
positive threshold d n . We consider an adequate kernel controlling the branching mechanism 
compatible with the discontinuous fragmentation kernel for the avalanches. We are able to 
follow the procedure stated in |BeDeLu 15]: first, we construct for each n a branching process 
on the finite configurations of E, greater than d n , having X n as base process (cf. Theorem 4.1; 
for a probabilistic description of this branching process see Remark 4.2), and then we project 
these processes on the set of all fragmentation sizes, in order to obtain the desired fragmentation 
process. We prove in Theorem 4.3 that this is actually a branching process with state space S^ 
and that starting from a sequence of fragmentation sizes (xk)k ^l £ S' , it is possible to restrict 
the process to the sequences of sizes of the form /3*(1 — fi^Xk, i, j G N, k ^ 1. This emphasizes 
that an avalanche has a fractal property, which is rather an infinite dimensional one, having its 
origin in the physical avalanche models studied in the papers mentioned in the first part of the 
Introduction; for more details see assertions (to) and ( iv ) of the Final Remark of the paper. 
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2 Fragmentation kernels for avalanches 

First order generators of jump processes. Let (E, dist) be a metric space, denote by B(E) 
the Borel cr-algebra of E, and by pB(E) (resp. bB(E), bpB(E)) the set of all real-valued positive 
(resp. bounded, positive and bounded) B(F , )-measurable functions on E. We endow bB{E) 
with the sup norm || • and denote by Ci(E) the set of all bounded Lipschitz continuous 
real-valued functions on E. 

Let N be a kernel on (. E , B(E)), JV^O. For all x E E we denote by N x the measure on E 
induced by the kernel N, 


N X (A) := N(1 a )(x) for all A E B(E). 


Suppose that 


( 2 , 1 ) 


dist(x, y)N x (dy) < oo for all x E E 


and define the first order integral operator N induced by N as 

Nf(x) := [ [f(y) - f(x)]N x (dy) V/ G bpB(E) \/x G E. 

J E 

The operator N is well-defined on the space Ci(E ): the function Nf is a B(E) -measurable 
real-valued function for every / G Ci(E). 

Recall that under various specific conditions (imposed on the kernel N, see [EthKu 86]) 
the linear operator N becomes the generator of a Markov process X = (X t )t^o with cadlag 
trajectories, such that N is the Levy kernel of X , i.e., 


E* X s )=E a 


S<t 



f(X s ,y)N Xs {dy)ds, 


o Je 


for all x G E and / G bpB(E x E), f — 0 on the diagonal of E x E. 

Examples. We present now several examples of such operators and processes, occurring in the 
study of the fragmentation phenomena and of the avalanches. 

1. The case of a continuous fragmentation kernel. Consider a fragmentation kernel 
F, that is, a symmetric function F : (0, l] 2 —> M + , and recall that F(x,y) represents the 
rate of fragmentation of a particle of size x + y into two particles of sizes x and y. If F is 
a continuous function then there exists an associate stochastic differential equation, as it is 
obtained in [FouGi 03] and [BeDeLu 15]. More precisely, assume that the fragmentation kernel 
F : (0, l] 2 —> M + is a continuous symmetric map. Moreover, F is supposed continuous from 
[0, l] 2 to M + U {+oo} and define the function 

- y{x - y)F(y,x - y)dy for x > 0, 
ip(x) = t x 


0 


for x — 0, 
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which is supposed continuous on [0,1]. The function ip(x) represents the rate of loss of mass 
of particles of mass x. With the notations from [FouGi 03] consider the operator T defined as 


Ff(x)= / [f{x-y) ~ f(x)]~ — -F(y,x-y)dy, x e [0,1]. 

Jo x 

We write J- as a first order integral operator. Let E = [0,1] and consider the kernel N F defined 
as 

(2.2) N F f(x):— [ f(z) — F(x — z,z)dz, x G E. 

Jo x 

The condition imposed to the function implies that N F satisfies the integrability condition 
(2.1) and one can see immediately that E = N F on CfE), that is, 

rf(x)= f\f(z)~ f(x)}N*(Az) V/ £ C'i(E) Vx£E. 

Jo 

Observe also that N F = 0 and 

(2.3) N F {dz) = — F(x — z, z)t(o jX )(z)dz, if x > 0. 

oc 

By Proposition 2.2 from |BeDeLu 15], using the existence and the uniqueness of the solution 
of the associated stochastic differential equation of fragmentation, under some mild conditions, 
N F becomes the generator of a Co-semigroup of contractions on C(E) and consequently of a 
standard (Markov) process with state space E. 

In general, we cannot associate such a stochastic differential equation of fragmentation be¬ 
cause we do not assume that F is continuous. However, in Section 3 below we succeed to write 
down a stochastic equation of fragmentation, appropriate to the avalanches. 

2. The case of a bounded kernel N. For the reader convenience we present now the 
classical situation of a bounded kernel (see, e.g., [EthKu 86], page 163), as we need to apply it 
to the fragmentations kernels for avalanches. 

Assume that N1 < oo and denote by A(x) the total mass of the (finite) measure N x , x G E, 

A(x) := Nl(x) G isl¬ 
and consider the induced normalized Markovian kernel N°, 

N° = -N. 

A 

Consequently we have 

(2.4) Nf(x) := AM f (/(a) - /M]iV°(da) V/ e bpB(E). 

J E 

Suppose that N is a bounded kernel. Then N becomes a bounded linear operator on bB(E), 

N = N — XI, 
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and it is the generator of a Co-semigroup (Pt)t^o on bB(E ): 

P t := e™, t ^ 0. 

Each P t is a Markovian kernel on E, more precisely, if we set 

(2.5) \ 0 : = ||iVl| joe and N' := + (1 - A)/, 

then P t f = e~ tXo ( A °j) N lk f, where N' k := N' o . . ,o N[ . With the above notations we 

k times 

have N = \ 0 N', that is 

Nf(x) = A of [f(y) - f(x)]N' x (dy) V/ G bpB(E) Wx G E. 

J E 

(2.6) The operator N is the generator of a (continuous time) jump Markov process (X t )t^o- 
More precisely, if v is a probability on E, we consider the time-homogeneous Markov chain 
{Y(k))k ^o in E, with initial distribution v and transition function N ', that is o T _1 (0) = v 
and 

W[f(Y(k + l))|y(0),... ,Y(k)\ = N'f(Y(k)) VfebpB(E). 

Let further (V t ) t ^\ be a Poisson process with parameter A 0 (i.e., P v (V t — k) — e~ tA ° , for 
all ken and t ^ 0) which is independent from (Y(k)) k ^o- Then the process (X t ) t ^ 0 defined as 

X t := Y(V t ), t jk 0, 

is a cadlag Markov process with state space E, transition function (Pt)t ^o and initial distribu¬ 
tion v. 


3. The case of (discontinuous) fragmentation kernels for avalanches. Consider again 
a fragmentation kernel F : (0, l] 2 —> R+. The following assumption is suggested by the so 
called rupture properties , emphasized in the deterministic modeling of the snow avalanches (cf. 
the papers mentioned in the Introduction). 


(Hi) There exists a function $ : (0, oo) —> (0, oo) such that 


F(x,y ) = <f> 



for all x, y > 0. 


Since the fragmentation kernel F is assumed to be a symmetric function, we have 





for all z > 0. 


Example of a fragmentation kernel satisfying (Hi). Fix a "ratio" r, 0 < r < 1, and 

consider the fragmentation kernel F r : [0, l] 2 —> R+, defined as 


F r (x,y) 


l(Mf) + Mf))M f x,y> 0, 
0 , if xy — 0. 
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We have F r (x,y ) = $ r (^) for all x, y > 0, where d* 7, : (0, oo) —» (0, oo) is defined as 

$ r (-) : = + S 1/r (z)), z> 0. 

Clearly, the function $ r is not continuous. Let ($())„ be a sequence of continuous functions 
such that ( < f>^-dx) ri is a sequence of probabilities on (0, oo), converging weakly to 5 r (dx). Let 


K(x.y) 



x, y > 0. 


The function F F is not symmetric, therefore we consider its symmetrization F£, 

K( x ,y) '■= 7 i (K(x,y) + K(y,x)), x,y> 0, 

and let N F ™ be the corresponding kernel given by (2.2). We have by (2.3) 

Nf”(dz) = —F^(x — z, z)t(o x -)(z)dz, x>0. 
x 

F r 

From the above considerations one can see that for each x > 0 the sequence of measures (N x n ) n 
converges weakly to ^-^-\fix5y x + (1 - j3)x8(\-p) x \, where /3 := 

Since the sequence of probabilities ( < fy/ r - dx)„ is approximating 5i/ r {dx) on (0, oo), it follows 

that ( N Fn/ ) n converges weakly to E-[/3x6p x + (1 — /3)x8(i-p) x ]. We conclude that the kernel 
N pr associated with F r is given by the following linear combination of Dirac measures: 


(2.7) N pr := \ 0 (/3x5p x + (1 - /3)x5 ( i-p) x ), 

where A 0 := 1 — 4 ^ . In this case the kernel N h is no more Markovian and has no density 
with respect to the Lebesgue measure. 


Further on in this paper we take 

E:=[ 0,1]. 

We fix now a sequence of thresholds for the fragmentation dimensions (d n ) n ^i C (0,1) 
strictly decreasing to zero. For each n ^ 1 let 

E n . [ d n , 1] and E n . d n ), Eq . E\. 

We assume that d\ < /3 and d n+ i/d n < f3 for all n ^ 1. 

Let n ^ 1 be fixed. Then E n = (Jfc=i F ' k _ v The kernel N pr from Example 3, given by (2.7), 
is used to define the kernel on E n as 

n 

Kf = 'E 1 K-, NFT (f 1 nJ Vf£tvB(E n ). 

k= 1 


Observe that, using (2.7), 

(2.8) the measure (N^) x is carried by [d k +i,x\ for every x E E' k , 
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and consider the corresponding Markovian kernel Nf, defined as in (2.5), 


1 / N r l\ 

Nnf = yKf+(l-^)f V/ g bpB(E n ). 

Following the procedure from Example 2, we may associate a first order integral operator 

r n f{x) := Nrf(x) = [ (f(y)-f(x)]N r n (dy) WfebpB(E n ) Vx e E n . 

J E n 

By (2.6) the operator T T n = A 0 (i\7 r — I) is the generator of a (continuous time) jump Markov 
process X n = (Xf) t>0 , defined as 

(2.9) X? := Y n (Vfi, t > 0. 

Its transition function is Pfi := e- 7 " 4 = e _A ° 4 'Ek^o n) k for ah t ^ 0 and let U n = (P„) a>0 

be the resolvent associated to the process X n , i.e., U™ = f^°e~ at Pfi d£, a > 0. Recall that 
(Y n (k))k ^o is the time-homogeneous Markov chain with transition function N'fi and (Vfit^o is a 
Poisson process with parameter A 0 . It is precisely the random walk described in Introduction. 


Remark 2.1. (i) Consider the integral operator P r defined as 

E r f(x) := N^f(x) = [ [f(y) - f(x)]N Fr (dy) f G bpB(E) x G E, 

that is, we take the kernel N pr instead of Nf. Then by the procedure described above there exists 
a jump Markov process X = (X t ) t ^ 0 , with state space E, having P r as generator, X t := Y(Vfi, 
t ^ 0, where (Y(k))k^o is the time-homogeneous Markov chain with transition function N’ r , 
where the Markovian kernel N ,r is obtained from N r as in (2.5). In the next section we show 
that the process X is related to a stochastic equation of fragmentation, associated to F r , the 
discontinuous fragmentation kernel for avalanches. 

( ii ) As in [BeDeLu 15], Proposition 5.2, we are able to solve the martingale problem asso¬ 
ciated to the bounded operator Pj (resp. P r ): for every f G bpB(E n ) (resp. f G bpB(E)) and 
each probability v on E n (resp. f G bpB(E)), the process 

f(X?) - fpjJiX^ds, (resp. f(X t ) - f T f(Xfids), t> 0, 

Jo Jo 

is a martingale under F u = JP x u(dx), with respect to the filtration of X n (resp. of X). In 
the Section j below X n (defined by (2.9)), will become the base process of a branching Markov 
process on the finite configurations of E n (see Theorem j.l below). 

For every iG [0,1] let 


E^ x -.= {nl-fiyx: i,j G N} U {0} 


and for n > 1 let 


J (3,x,n 


Ep, x n E n 
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Proposition 2.2. The following assertions hold for the conservative Markov process X n with 
state space E n , n^l. 

(i) Ift^s then a.s. X™ ^ Xf. 

(ii) The sets [d n ,x] and [d n ,x) are absorbing subsets of E n for every x G E n with respect to 
the resolvent U n of X n . In particular, Ef^ is an absorbing subset of E n . 

(Hi) If x G E n then P x -a.s. Xf G Ep xn for allt^O. 

H If(P t n ) tj.o is the transition function of X n , then 

Pf) x l = PJf x for all and x G E n . 

Proof. ( i ) Assume that s = 0 and let x G E n . Then Xq = x P x -a.s., so, we have to prove that 
P ,x (Xf G [d n ,x]) = 1, or equivalently, that Pff( l( X) i])(x) = 0 for all t. Observe hrst that if / = 0 
on [d n ,x] then by (2.8), Nfff = 0 on [0, x] and by induction we get (N'ff) k f = 0 on [d n ,x\ for 

all k^O. Therefore Pff = e~ Xot Ylk^o ^°k\ (N'n) k f a l so vanishes on [d n ,x\ for all f^O. The 
case s > 0 follows using the Markov property of X n . 

(ii) We argue as in the proof of Proposition 3.1 from |BeDeLu 15]. By the above consider¬ 
ations we have P x -a.s. Xf ^ x for all t ^ 0 and therefore the entry time D( Xi i] of the set (x, 1] 
is inhnite. Indeed, we have hrst D/ y ^ = oo P y -a.s. for all y G E n and therefore, if d n ^ y < x, 
then [-D( yi i] = oo] C [-D(x,i] = oo], P y ([H( Xil ]) = oo]) = 1. By (A.1.2) from [BeDeLu 15] we 
conclude that [d n ,x] is absorbing. The set [0, x) is also absorbing as a union of a sequence of 
absorbing sets, [d n , x) = (j k [d n , d n + (x - d n - |)+]. 

(Hi) It follows by induction from (2.7) that the probability measure (N!f) k is carried by 
Ep tXtU for every x G E n and k ^ 0. Hence P x (Xf G Ep jXjn ) = Pf l (l Ef , xn )(x) = 1. 

(iv) It is sufficient to show that for every / G B(E n+1 ) and k ^ 1 we have on E n (Nf +l ) k f = 
(N!f) k (ft En ). Indeed, one can see that N x +1 f = N x + 1 (ft En ) = Nf(ft En ) on E n and therefore, 
by induction, (N p +1 ) k f = (K +l ) k (fl En ) = (K) k (fl En ) on E n . □ 

If X n , n^l, is the Markov process with state space E n , given by (2.9), and having (Pf 1 )^ o 
as transition function, constructed from the discontinuous fragmentation kernel F r , then by 
Proposition 2.2, assertions (ii) and (Hi), conditions (H-f) and (H 3 ) from [BeDeLu 15], Section 
3, are fulfilled by the processes X n . 

3 The corresponding stochastic equation of fragmentation 

To emphasize the stochastic equation of fragmentation which is related to our stochastic model 
for the avalanches, we rather consider for simplicity the kernel N pr on E instead of N x , n ^ 1, 
and the associated pure jump process X = (X t ) t ^ 0 with state space E, given in assertion (i) of 
Remark 2.1. 


We state now the stochastic equation of fragmentation for avalanches: 


X 2 a _da 


(3.1) 


X t = X 0 — 2X 0 j3(l — j3) 


0 




where p(da, du, ds) is a Poisson measure with intensity q := dadwds. 

The next theorem gives the claimed relation between the equation (3.1) and the process X. 

Theorem 3.1. The stochastic equation of fragmentation for avalanches (3.1), with the initial 
distribution S x , x G E, has a weak solution which is equal in distribution with (X, P x ). 

Proof. Define the bounded kernel K on M by 



Using (2.7) one can see that 



where on left hand side the kernel N pr occurring in the definition of J- r is extended from E to 
M with zero on the complement of E. In particular, we have Kf(x) = / f(y — x)N pr (dy) for 


x G [0,1]. By a strait-forward procedure (see e.g., Lemma (2.2) of Chapter V from [BIGe 68]), 
applied to the measure N pr which is carried by (0, x), we get for all i,m,s6 M: 



(3.3) 


where 


inf{u > 0 : N pr ((0, u]) > s} — x if 0O,0 ^m^1,i 6 [0,1], 
0 else, 



with the convention /(oo) = 0. 

Let further (/C, C£ (M)) be the operator defined as: 


/C/(a;) = b(x)f'(x) + [ [f(x + y)~ f(x) - yt { \ y \< l} f(x)}K x {dy), 

Jr 


cf. [Jac 79], page 434. Taking b(x) = J R yK x (dy) = 2A 0 f3(/3 — l)x 2 if x G [0,1] and b(x) = 0 
elsewhere, we clearly have T r = /C. 

We consider the following stochastic differential equation, applying a method from [Jac 79], 
page 479: 


AX, = 2X„P(P - l)A' f 2 _df + w t dp t - .-idq,. 


(3.4) 


where w t is a process that depends on X t , such that Kf(X t _ ) = J p f R dudsf(w t (-, u, s))l[u) t (.,u,s)^o]- 


Note that in general, the existence of the process Wt follows by Theorem (14.53) from [Jac 79], 
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applied to the measure duds. However, the main observation here is that using (3.3), we can 
take w t (-, u, s ) = r(X t _,u, s ), or explicitly 

W t (-,U,S ) = l[0<,t<i](u)l[0<X t _CL] ((/3 - l)X t _l[^_ <Xt _] - px t _ + OOl[ A - t _^]) • 

From the above considerations, the stochastic differential equation (3.4) may be rewritten in 
the form (3.1). 

By assertion (ii) of Remark 2.1 we know that the martingale problem associated to the 
bounded operator T r has a solution. Consequently, Theorem (14.80) from [Jac 79], page 48, 
implies the existence of the solution of the stochastic differential equation (3.4). □ 

4 Branching and fragmentation processes of an avalanche 

For a Borel subset E of [0,1] define the space of finite configurations of E which is the following 
set E of finite positive measures on E: 

E 5 Xk : ko G W,x k G E for all 1^/c^/co > U {0}, 

Ifcyfco J 


where 0 denotes the zero measure. We identify E with the union of all symmetric m-th powers 
E ( ' m ' ) of E: E = (Jm^o E^ m \ where E^> := {0}; see, e.g., [INW 68], [BeOp 11], [BeLu 15], and 
[Si 68]. The set E is endowed with the topology of disjoint union of topological spaces and the 
corresponding Borel a- algebra 13(E). 

If p\,p 2 a re two finite measures on E , then their convolution p\ *p 2 is the finite measure on 
E defined for every h G pB(E ) by 


pi * p 2 (du)h(u) ■= / pi(dz/i) / P2(dv 2 )h(vi + v 2 ). 


If ip G pB(E), define the multiplicative function Cp : E — y M + as 

{ \v{xk), if x = (x k ) k>1 Gh,x^ 0, 

k 

1 , if x = 0. 

Recall that a bounded kernel N on E is called branching kernel if 

N, 1+1J = N^* N u for all p,v G E, 

where Ay denotes the measure on E such that f gdN^ = Ng(p) for all g G pB(E). Note that 
if N is a branching kernel on E then A r 0 = <5 0 G M(E ). A right (Markov) process with state 
space E is called branching process provided that its transition function is formed by branching 
kernels. The probabilistic description of a branching process is the following: if we take two 
independent versions X and X' of the process, starting respectively from two measures p, and 
//, then X + X' and the process starting from p + p! are equal in distribution. 
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Branching processes on the finite configurations of E n . For all n ^ 1 we define the 
Markovian kernel B n from E n to E n as 

(4.1) B n h(x ) := —Y lE Li( x ) d k h (y, y)y( x - y), h e bpB(E n ), x e £ n , 

Ci ( ) _ 

v 7 E/3 jX 3y<x 

where for d > d n and g £ 1]) we consider the function £ bpB(E n ), the extension of g 

to E n with the value g(d) on [ d n , d), 

dg{y) ■= g{d)Md n ,d)(y) + 9(y)Md,n(y)^ y e E n, 

and a(x) := Y^ y(x — y) < oo for all a: £ E n . 

Ep lX ^y<x 

Observe that sup^g^ B n li[x) = 2, where for a function / £ pB(E n ) we consider the mapping 
If : E n —>■ M + defined as lf(p) '■= f f dp, p £ E n . Therefore we my apply Proposition 4.1 from 
|BeLu 15], to construct a transition function (Pp)t ^o on E n , induced by {P™)t^ o> the transition 
function of the Markov process X n on E n , given by (2.9), and by the kernel B n . The existence 
of a branching process with state space E n and transition function (PJ 1 )^o is given by the first 
assertion of the following theorem. It is a version of Proposition 5.1 from [BeDeLu 15]. 

If X \ 1 ... , Xk £ E and x = S X1 + . .. + S Xk £ E, we put 


k 



If n ^ 1 then let 

k 

E 0,x,n ■ Ef),Xj,n- 

3 =1 

Theorem 4.1. Let n ^ 1. Then the following assertions hold. 

(i) There exists a branching standard (Markov) process X n = (Xf) t f 0 . induced by (P™)t^ o 
and the kernel B n , with state space E n and having the transition function (Pf 1 )^ o- 

(ii) For every x £ E n , y £ Ep^, n , and t ^ 0 we have P y -a.s. Xf £ Ep iXjU . 

Proof, (i) We argue as in the proof of Proposition 5.1 from [BeDeLu 15]. Consider the vector 
space C n defined as 


C n ■— {/ : [d n , 1] —> M : /|b' £ C(E' k ) such that lim f{y) £ R for all k — 0,..., n — 1}. 

k y/'dk 

We claim that (P™)t^o induces a Co-semigroup of contractions on C n . Indeed, we observe 
first that the Markovian kernel N(f becomes a bounded contraction operator on C(E n ) and 
consequently P f n (C n ) C C n for all t^O. Also, for all / £ C n we have \im t \o Pff — f uniformly 
on E’ k for each k — 0,..., n -< 1. Note that the kernels B n have the following property 

if ip £ pB(E n ), then B n (p £ C n . 

Now condition (*) from the proof of Proposition 5.1 from |BeDeLu 15] is verified by the vector 
space C n . So, there exists a standard process with state space E n , having (Pf 1 )^o as transition 
function. 
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(ii) Let ip := II Ep x n and observe that (p = l^r-— ■ Since for y G Ep xn we have Ep t y>n C 
Ep iX>n , we only have to prove that P[ l (p(x) = 1. By Proposition 4.1 from |BeLu 15] we have 
Pfip = h t (ip), where h t ((p) h t G bpB(E n ) is the unique solution of the integral equation on 
E n ^ 

h t (x) = e~ t P[ l ip(x) + f e~^~ u ^ P[f u B n h u (x)du, t ^ 0, E n . 

Jo 

It is sufficient to show that 

(4.2) = 1 on E 0 ,x,n for every x G E n . 

Indeed, if x G E n , x = 8 X1 + ... + S Xk , since the function <p i— y h t (ip)(x) is increasing, we have 

k k 

P/V(x) = J ~[ht((p)(xj) >'Wh t il B p^ n )ixj) = 1, 
i = 1 i =1 

where the last equality holds by (4.2). To prove it let x G E n and assume further that ip := 
1 Ep xn - Recall that (cf. the proof of Proposition 4.1 from [BeLu 15]) ht is the point-wise limit 
of the sequence (/i™) m eN defined inductively as follows: = e t PJ L ip + /Je ^ ^ P™_ u B n <pdu, 

K +l = e~ t P^ l ip + fo e-( t ~ u )P£_ u B n h™du, m ^ 0. Assertion (in) of Proposition 2.2 implies that 
PJip = 1 on Ep x n for all t ^ 0 and by (4.1) we also have B n (p = 1 on Ep jXjn . It follows that 
on Ep, xn we have h® = 1 and by induction h™ = 1 for every m ^ 0. We conclude that (4.2) 
holds. □ 

Remark 4.2. (i) The branching process constructed in Theorem 4-1 on E n has the following 
description: An initial particle starts at a point of E n and moves according to the base process 
X n , until a random time, when it splits randomly into two new particles, its direct descendants, 
placed in E n . Each direct descendant starts from a position dictated by the "branching kernel" 
B n and moves on according to two independent copies of X n and so on. 

(ii) In the above description the number of direct descendants is two since the kernel B n 

-—— /qN 

from E n to E n , defined by (f.l), is actually carried by En . Accordingly, a binary fragmentation 
process will be involved. 

(Hi) In contrast with the situation from [Si 68 ], the starting position of a descendant may 
be different from the terminal position of the parent particle. This property is compatible with 
our aim to count with the branching processes the number of fragments of the avalanche, having 
the same size. So, since for our model a particle represents a size, then clearly the descendants, 
resulting after a splitting, should represent sizes which are sm.aller than (in particular, different 
from) the parent size. 

Branching processes on the space of fragmentation sizes. As in |BeDeLu 15], we intend 
to construct a process with state space the set S 4 of all decreasing numerical sequences bounded 
above from 1 and with limit 0, 

S l := {x = (x k )k >l Q [0,1] : (x k )k>i decreasing, lim x k = 0}. 

k 

Recall that a sequence x from S^ may be considered as "the sizes of the fragments resulting 
from the split of some block with unit size" (cf. |Ber 06], page 16). It is convenient to identify 
a sequence x = (x k )k^i from S' 4 with the cr-finite measure /i x on [0,1], defined as 
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j > if x ^ °> 

hx • \ k 

[ 0 , if x = 0, 

where the zero constant sequence 0 is identified with the zero measure, /io = 0. 

Let further 

S := {x = G 5''*' : 3 G N* such that Xk 0 > 0 and Xk = 0 for all k > k 0 }. 

The mapping x i—» p x identifies S with lj n>1 E n \ {0}. For x 6 we write x = p x where 
it is necessary to emphasize the identification of the sequence x with the measure /i x . 

In order to consider branching kernels on S^ we need to have a convolution operation between 
finite measures on S'h We first endow with a semigroup structure: if x, y G then the 
sequence x+y G is by dehnition the decreasing rearrangement of the terms of the sequences x 
and y. The convolution may be now introduced as in the case of the space finite configurations: 
if pi,P 2 are two finite measures on S^, then their convolution p\ * p 2 is the finite measure on 
S l , defined for every h G pB(S l ) by: J sl px * p 2 (du)h(u) := f si p,(dzq) J sl p 2 (du 2 )h(u 1 + v 2 ). 
The branching kernels on and the branching process with state space S^ are now dehned 
analogously, and the probabilistic interpretation remains valid. 

Define the mapping a n : S^ i —> E n as a n (x) := p x \ b„, x = p x G S^. We have a n (0) = 0 
and a n Igs = Idgr. Dehne also 

1 Ej-n h/n 

Soc ■= {(x n )n^i G E n : x n = a n (x m ) for all m > ri^l}. 

n^l 

By Proposition 4.5 from |BeDeLu 15] we have: 

(4.3) The mapping i : S^ —> Soo, dehned as «(x) := (a n (x)) n ^i, x G S^, is a bijection. 

Using assertion (• iv ) of Proposition 2.2, we may apply Proposition 4.6 from [BeDeLu 15] to 
get the following result: 

(4.4) Let x G and x n := cc n (x) G E n ,n^l. If t > 0 then the sequence of probability 

measures (Pt,x n )n ^l is projective with respect to (E n ,a n ) n ^i, that is P’" Xn 1 +1 ° ot~ l = P" Xn for 
all n^l. 

By the identihcation of with S 0Q , given in (4.3), the projective system of probabilities 
from (4.4) may be used to apply Bochner-Kolmogorov Theorem (see, e.g., [BeCi 14]) in order 
to induce a transition function on the space of all fragmentation sizes, as in Proposition 4.7 
from [BeDeLu 15]: 

(4.5) There exists a Markovian transition function (Pt)t^o on such that for each x G and 
n^l we have P t , x ° a” 1 = P” Xn , where x n := a n (x). 

We can state now the result on the fragmentation-branching processes related to the avalanches, 
having as state space the set 5'"*' of al fragmentation sizes, endowed with the topology induced 
by the identihcation from (4.3) with S oo (equipped with the product topology). 

Theorem 4.3. The following assertions hold. 
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(i) There exists a branching right (Markov) process X = (X t ) t ^ 0 with state space S^, having 
cadlag trajectories and the transition function (Pt)t^o, given by (f.5). 

(ii) If x G [0,1], y G S^, y^x, and tjz 0 we have X t ^x P y -a.s., where if y = (yk)k^i, y^x 
means that yk^x for all k^l. 

(Hi) For each xe£, the set S^ x := {y = (pk)k^i G : yk G Ep tX for all 1} is absorbing 
in S^, that is, ify G Sp x then P y -a.s. X t G Sp x for all t ^ 0. 

Proof. We prove that each kernel P t , t ^ 0, is a branching kernel on S'h Let x, y G S^, we have 
to show that 

Pt,x+ y (g) = Pt, x * Pt, y (g ) for a11 9 e bpB(S l ). 

Since by (5.10) from [BeDeLn 15] we have B(S^) = a (\J n ^{/ ° a n '■ f € bpB(E n )} S j, we may 

assume that g = f o a n for some n ^ 1, with / G bpB(E n ) and using (4.5) it remains to prove 
that 


(4.6) 


pn 

1 £,x n +y n 


(/) = Pt 


t,x„ 


* p t?y„(/) for a]1 / e bpB(E n ), 


where x„ := a n (x). Here we used the equality a n (x + y) = a n (x) + a n (y), where the second 
sum is the usual addition of measures from E n . But the equality (4.6) is precisely the branching 
property of the kernel Pf on E n , which clearly holds by Theorem 4.1. 

The existence of the process X claimed in assertion (i) and assertion (ii) are consequences 
of |BeDeLu 15], Theorem 5.3 and Corollary 5.5 respectively. Note that a key argument in the 
proof of the path regularity is the existence of a compact Lyapunov function (a superharmonic 
function having compact level sets); see [BeLuOp 12] and also |BeRo 11] for some related 
results. 

(Hi) By the right continuity of the process X and since the set Sp x is closed, it is sufficient 
to show that for each t ^ 0 we have P t ( )(y) = 1 for every y G S^ x . We have: y G Sp x if 

and only if a n (y) G Ep^, n for all n ^ 1. On the other hand by assertion (ii) of Theorem 4.1 
we get Pf(ts —- )(y n ) = 1 provided that y n belongs to FL xn . Using again (4.5), we conclude 

■Pfl ,x.,n ' 

that 


)(y) = P t , y(p| = lim Pt,y K 1 (Ep, x ,n) ) = lim P t n yn (Ep^ n ) = 1, 


n> 1 


where y n = a n (y) G Ep^ n . □ 

Final remark, (i) Recall that at the origin of the equation (3.1), the stochastic equation 
of fragmentation for avalanches, is the discontinuous fragmentation kernel F r . This equation 
should be compared with the stochastic differential equations of fragmentation with continuous 
fragmentation kernels from |FouGi 03] and equation (2.4) from |BeDeLu 15]. 

(ii) The statement of assertion (i) of Theorem 4.3 is valid, with the same proof, for the 
fragmentation process constructed starting with a fragmentation equation (with a continuous 
fragmentation kernel). More precisely, the right (Markov) process, given by Theorem 5.3 from 
|BeDeLu 15], is in addition a branching process with state space This is similar to the 
branching property proved for the fragmentation chains in Proposition 2.1 from [Ber 06]. 
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(Hi) Since by assertion (Hi) of Theorem 4.3 the set Si x is absorbing, it is possible to restrict 

the fragmentation-branching process X to this set; for the restriction procedure see, e.g., (12.30) 
in [Sh 88] and Appendix A.l in |BeDeLu 15]. 

(iv) Assertion (Hi) of Theorem 4.3 emphasizes the fractal property of an avalanche, claimed 
in the Introduction and closed to its real physical properties: if we regard the fragmentation¬ 
branching process on the set Sp x (which is possible by restriction, according to the previous 
assertion (Hi)), then independent to the sequence of sizes x of the initial fragments, from the 
moment when the avalanche started, and remaining constant in time, the ratio between the 
resulting fragments are all powers of (3. 

(v) In a future paper we intend to give a probabilistic numerical approach to the avalanches, 
based on the results from this paper, and using appropriate stochastic equations of fragmenta¬ 
tion. 
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